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SUMMARY 


A computer program written to calculate the static pressure position error 
of airspeed systems contains five separate methods for determining position error, 
of which the user may select from one to five at a time. The program uses data 
from both the test aircraft and the ground-based radar to calculate the error. In 
addition, some of the methods require rawinsonde data or an atmospheric analysis, 
or both. 

The program output lists the corrections to Mach number, altitude, and static 
pressure that are caused by position error. Reference values such as angle of 
attack, angle of sideslip, indicated Mach number, indicated pressure altitude, 
stagnation pressure, and total temperature are also listed. 


INTRODUCTION 


The measurement of true airspeed and altitude is extremely important when 
flight testing aircraft. To increase the accuracy of these measurements, the 
static pressure position error must be determined, then accounted for. Though 
many methods have been developed and used successfully, the data analysis can be 
tedious. To decrease the workload, a computer program (PERROR) was written that 
contains several methods (all incorporating radar measurements) for determining 
static pressure position error for a given air data system. The program accesses 
data from a file that contains ground-based radar and test aircraft information. 
Additionally, atmospheric analysis data are entered in table form if available. 
PERROR then calculates the position error for the system. 

This report describes the methods that are used in the program, the inputs 
necessary for each method, and the limitations of the methods. In addition, the 
output from an actual test is presented. 


ALPHANB 

BETANB 

CP 

DANG 

DESCP 

DESCT 

DHP 

DHPG 


NOMENCLATURE 

angle of attack, deg 

angle of sideslip, deg 

pressure coefficient, (PICT - P)/QCIT 

difference between azimuth angle and direction of pressure gradient 
descent pressure method 
descent temperature method 

calculated difference between true and indicated pressure altitude 

difference between geometric and pressure altitude as a function of 
horizontal position 


DM 


calculated difference between true and indicated Mach number 


DPR 

DR 

OZ 

DZEN 

DZES 

DZH 

DZHTABL 

E 

FLIGHT 

G 

GGHTBL 

GH 

HP 

HPREK 

HPT 

h 

IETSV 

IET 

IETAD 

II 

1ST 

ISTAD 


calculated difference between true and indicated static pressure 
distance from radar site to aircraft, ft 

constant adjustment to radar altitude (ft) that allows the computed 
pressure altitude at the reference point to equal the calibrated 
pressure altitude; used if KK * 0 and LL * 0 

difference between geometric altitude and indicated pressure altitude 
as a function of elevation angle along a survey run 

difference between geometric altitude and indicated pressure altitude 
at the maximum elevation angle along a survey run 

difference between geometric and pressure altitude as a function of 
geometric altitude 

table of pairs of DZH as a function of altitude; used if KK * 0 and 
LL £ 0 

radar antenna elevation angle, deg 
flight number 

pressure gradient obtained from atmospheric analysis, ft/n.mi. 

table of triples of G and GH as a function of altitude; used if KK * 0 
and ISURVEY = 0 

direction of the pressure gradient (direction of decreasing Z - HP at 
constant geometric altitude, measured in degrees from true north) 

* 

true pressure altitude, ft 

reference pressure altitude, ft; used if MM # 0 
test or indicated pressure altitude, ft 
geopotential altitude, ft 

end time (HR, MIN, SEC, MSEC) of the survey run; used if ISURVEY t 0 
stop time (HR, MIN, SEC, MSEC) 

end time (HR, MIN, SEC, MSEC) of the acceleration-deceleration run 
associated with a survey run; used if ISURVEY * 0 

flag set to nonzero to calculate position error using radar-rawinsonde 
method 

start time (HR, MIN, SEC, MSEC) 

start time (HR, MIN, SEC, MSEC) of the acceleration-deceleration run 
associated with a survey run; used if ISURVEY * 0 
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ISTSV start time (HR, MIN, SEC, MSEC) of the survey run; used if ISURVEY * 0 

I SURVEY set to nonzero for a survey run (a survey run is used only when KK * 0) 


I UNITS 


KK 


TiEVD 


LL 


M 

MI 

MM 


NDZH 


NGGH 


NN 


OO 


set to 0 indicates standard English units for all inputs; set to 1 

indicates pressure units are in lb/in 2 and other units are in stand- 
ard English units; set to 2 indicates pressure is in N/m 2 and TT is 
in kelvins 

flag set to nonzero to calculate position error using level 
acceleration-deceleration method 


level acceleration-deceleration method 

flag set to nonzero to calculate position error using descent pressure 
method 

Mach number 

indicated Mach number 


£ lag set to nonzero cu odiculdte posi ti on error 
perature method 


using descent tern- 


two times the number of pairs of entries in the DZH table; used if 
KK * 0 or LL * 0 


three times the number of triples of entries in the GGH table; used if 
KK * 0 and ISURVEY = 0 


flag set to nonzero to calculate position error using total temperature 
method 


flag set to 1 to compute PT from QCIT and PICT from the merged file; 
set to 2 to compute PT from VCT and HPT from the merged file; default 
of 0 will use PT as it comes from the merged file 


ambient stati n pressure , lb/f t 2 


PICT indicated static pressure, lb/ft 2 

PR ambient pressure determined from rawinsonde balloon measurements, lb/ft 2 

PT indicated total pressure, lb/ft 2 


QCIT indicated impact pressure, lb/ft 2 

qq flag set to nonzero to output pressure coefficient CP; default of 0 

will output pressure altitude correction DHP 

rra radar parameter identification for azimuth angle, deg 
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RRH 

RRTA 

RRX 

RRY 

RUN 

T 

TOT 

TS 

TT 

VCT 

Z 

Y 

Suffixes : 

DP 

I)T 

LD 

R 

TT 


radar file parameter identification for geometric altitude, ft above 
sea level 

radar file parameter identification for ambient temperature, °F 

radar file parameter identification for x coordinate, centered at sta- 
tion, positive east, ft 

radar file parameter identification for y coordinate, centered at sta- 
tion, positive north, ft 

run number 

ambient temperature, °F 

total temperature method 

standard temperature, °F 

total or stagnation temperature, °F 

calibrated airspeed, knots 

geometric altitude, ft 

specific heat ratio 


output calculated using 
output calculated using 
output calculated using 
output calculated using 
output calculated using 


the descent pressure method 

the descent temperature method 

the level acceleration-deceleration method 

the radar-rawinsonde method 

the total temperature method 


DESCRIPTION OF METHODS 


Level Acceleration-Deceleration Method 

The level acceleration-deceleration method involves deriving true pressure 
altitude HP from geometric altitude Z and a calculated or measured difference in 
geometic altitude and pressure altitude, Z - HP. Ideally, the test aircraft flies 
a constant-altitude acceleration-deceleration run to cover the desired Mach range. 
However, if desired, this method can be used for a climbing or descending maneuver. 
Geometric altitude is obtained from ground-based radar. 
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An atmospheric analysis is conducted to determine values of Z - HP, pressure 
altitude gradient G, and ambient temperature T. When the value of Z - HP is 
subtracted from the radar altitude, true pressure altitude is obtained. The value 
Z - HP is referenced vertically, directly above the radar site. The pressure alti- 
tude gradient is the slope of the pressure altitude in the horizontal plane and is 
represented in vector form. The magnitude is in geopotential feet per nautical mile 
of horizontal distance and in the direction of decreasing Z - HP. The origin of 
this vector is at the radar site. 

The atmospheric analysis is based on information received from rawinsonde 
soundings. NASA Ames Research Center, Dryden Flight Research Facility (Ames- 
Dryden) receives information from the Air Weather Service at the Air Force Flight 
Test Center, Edwards AFB. The Air Weather Service gathers data from rawinsonde 
soundings made twice a day at numerous locations. The analysis is made from 
constant-pressure-level maps from which pressure-height, temperature, and wind 
information can be extracted for the flight test runs. In addition, an Edwards 
rawinsonde sounding is generally scheduled by Ames-Dryden to provide data in the 
vicinity of the calibration flight, at approximately the time of the flight. Data 
from this sounding provide more points per unit of altitude than does the other 
sounding information. 

When using the level acceleration-deceleration method, two options are 
available to compensate for the horizontal pressure gradient: the nonsurvey 

option, using an atmospheric analysis, and the survey option, using a constant- 
speed survey run made along the same path as the acceleration-deceleration. The 
option descriptions follow. 

Nonsurvey option . — The nonsurvey option of the level acceleration-deceleration 
method uses both Z - HP and pressure-altitude gradient from the weather analysis. 

At each time interval (typically every second), a value of Z - HP as a function of 
vertical height DZH and a value of Z - HP as a function of horizontal position DHPG 
are calculated. The value of DZH is obtained directly from the weather analysis ; 
DHPG is calculated by multiplying the distance to the aircraft from the radar site, 
DR, by the pressure gradient G and the cosine of the difference between the azimuth 
angle and direction of the pressure gradient, DANG: 

DHPG = (DR) (G COS DANG) 

A schematic of this method is shown in figure 1(a). True pressure altitude is then 
determined by the equation 


HP = Z - DZH + DHPG - DZ 

where DZ is an adjustment. The variable DZ defaults to zero if there is no 
existing static-pressure position error correction. If a position error is firmly 
established for a specific flight condition, DZ can be used so that the calculated 
pressure altitude at the same flight condition is equal to the known pressure alti- 
tude. The previous equation can also be used for the "bootstrap" method where, 
again, the static pressure position error is known at a reference point, but values 
of DZH and DHPG are not used. In such a case, DZ is evaluated at the reference 
point. Finally, in all cases, DZ is applied as a constant throughout the entire 
test run. 
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Survey option . — With the survey option, DZH is used to provide a vertical 
component of Z - HP, as in the nonsurvey method* However, a different scheme is 
used to arrive at a horizontal component* During the flight test, a constant- 
Mach run is made at the same altitude and over the same path as the acceleration- 
deceleration run* From the survey run, a correction for pressure gradient error 
and radar error (such as, in elevation angle) can be obtained as a function 
of radar-antenna elevation angle. An initial value of Z - HPT = DZES is deter- 
mined at the highest elevation angle, where radar errors are at a minimum. A table 
of Z - HPT values (DZEN) as a function of elevation angle is generated from the 
survey run* Then a correction for each point along the acceleration-deceleration 
can be determined from DZEN at the corresponding elevation angle. This method is 
shown schematically in figure 1(b). The following equation is used to calculate 
pressure altitude: 


HP = Z - DZH - (DZEN - DZES) - DZ 


Radar-Rawinsonde Method 

The radar-rawinsonde method determines true Mach number directly. Adiabatic 
equations relating Mach number to stagnation and static pressure are used: 

M = [5(PT/P)°* 2857 - 5]°* 5 M < 1.0 

M = [(1.42851 - 0. 3571 43y - 0.0625y2 - 0.025y3 - 0.012617y4 

- 0.007l5y5 - 0.004358y 8 - 0.0087725y 9 )/y]0* 5 M > 1.0 


where 


y = 1 .839371/(PT/P) 

The supersonic equation accounts for the total pressure loss across a normal shock 
wave. The approximation is given in appendix A. It is assumed in the equations 
that there is no position error in the measured stagnation pressure. The static 
pressure as a function of geometric altitude is obtained from the rawinsonde 
listing for the flight test. 


Descent Pressure Method 

The descent pressure method, which can also be used for climbs, is similar to 
the level acceleration-deceleration method in that it calculates true pressure 
altitude by subtracting Z - HP from radar altitude. However, with this method the 
test aircraft is flown in a descent or ascent, as opposed to level flight, and no 
corrections are made for horizontal variations in the atmosphere. 

As in other methods, Z - HP is obtained from an atmospheric analysis, and DZ 
can be used to obtain the correct pressure altitude at a reference point in the run 
for which the static-pressure position error is known. As previously described, 
this correction is then applied throughout the rest of the run. 

HP = Z - DZH - DZ 
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Descent Temperature Method 


The descent temperature method, which can also be used for climbs, calculates 
true pressure altitude from stagnation temperature, indicated Mach number, and 
radar altitude. The difference between true and indicated pressure altitude then 
determines the correction for static-pressure position error. 

Procedurally, the calculation begins at a reference point (either the first or 
last point of the run) for which the true pressure altitude HPREF is required. 

The departure of the radar altitude from the reference altitude is divided into 
increments, each being converted to a pressure-altitude increment by the following 
method. First, the geometric altitude increments are converted to geopotential 
altitude by using the approximation 

Ah - (20,825,000 AZ) / ( 20 , 825 , 000 + Z) ft 

The resulting geopotential increment is converted to the associated pressure- 
altitude increment by the equation 

A HP - TS/T Ah 

In this equation the temperatures are calculated by averaging the temperatures at 
the end points of the altitude increments. The standard temperature TS is calcu- 
lated as a function of geopotential altitude. The ambient temperature T is deter- 
mined from the approximation 

T = TT/(1 + 0.2M 2 ) 

Iteration is required in this equation since Mach number is used. Indicated 
Mach number is first used to obtain T, from which position error (and hence a new 
Mach number) is calculated. This new Mach number is used in the preceding equa- 
tion for another iterative cycle. Convergence on true Mach number occurs rapidly. 

After the geometric altitude increments have been converted to pressure- 
altitude increments, they are summed to determine position error correction for the 
test point: 

y 

HP = HPREF + £ ( AHP ) j - HPT 

j=0 

where y is the number of altitude increments and the sum of the first two terms on 
the right-hand side of the equation is the true pressure altitude. 

Total Temperature Method 

The total temperature method solves for Mach number from the equation 

TT = T(1 + 0.2M 2 ) 


j 
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The specific heat ratio is assumed to be 1.4. The stagnation temperature TT is 
measured onboard the test aircraft, and the ambient temperature T is obtained from 
the atmospheric analysis, or directly from the rawinsonde. 

The correction for Mach number determined by this method is for total pitot- 
static position error, not just the correction for the static-pressure posi- 
tion error. 


METHOD LIMITATIONS 


General Limitations 

This program has several limitations. No data smoothing is done nor does the 
program seek and discard wild points. There are no provisions to compensate for 
pneumatic lag of the airspeed system. In addition, there are no corrections to the 
position error for real gas effects. One safeguard, however, is that if the value 
PT < P, calculations for that point are terminated. 

Radar Data 

The radar data used by Ames-Dryden are obtained from an FPS-16 radar system. 
Following a series of flight experiments, it was determined that the radar data may 
be questionable when the radar antenna elevation angle is below 7°. This is due in 
part to the large uncertainty of the refraction correction that exists at the low 
elevation angles. However, radar data from the lower elevation angles may be used 
if the survey option is exercised because the option will account for systematic 
radar errors. In addition, loss of beacon track will give unusable radar data. 

Level Acceleration-Deceleration 

When using the level acceleration-deceleration method with the nonsurvey 
option, the results are very dependent on the accuracy of the atmospheric analysis. 

Temperature Descent Method 

when exercising the temperature descent option, the altitude intervals should 
be less than 100 ft to adequately approximate the ambient temperature. 

PROGRAM USE 


General Description 

To implement the program the user must first create a merged file of the neces- 
sary aircraft and radar data. All other necessary inputs are contained on a single 
data card that must be submitted with PERROR. The program then reads the data for 
the specified time interval, and unit conversions of inputs to the standard English 
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system are made if required. If total pressure PT is not available directly, it 
is calculated from either impact pressure QCIT and static pressure PICT or veloc- 
ity VCT and pressure altitude HPT. The output format for the flight conditions 
are then written. 

At this point the corrections for position error are calculated, using any or 
all of the five methods, but usually no more than three. Each method uses common 
subroutines for calculating Mach number, pressure altitude (1962 Standard 
Atmosphere), and static pressure, for example. The main program calculates and 
writes the corrections for Mach number DM, static pressure DPR, and pressure alti- 
tude DHP caused by position error for each method used. After the calculations are 
performed for one time point using all methods selected, the program completes the 
calculations for the next time point. A listing of the program is presented in 
appendix B. 


Program Inputs 

A block diagram of the input sources to PERROR is shown in figure 2. In addi- 
tion, table 1 lists the necessary input parameters for a particular method. First, 
a flight data base file and a radar data file are created. Rawinsonde data may be 
included when creating the radar data file* The data card required to run PERROR 
must contain the following quantities: FLIGHT, RUN, IUNITS, OO, and QQ. Zero is 
the default value for the last four quantities* The indices II, KK, LL, MM, and NN 
are used to designate which position error methods are to be used (refer to the 
NOMENCLATURE section). 

The following is an example of a data card used when running PERROR (the DZH 
and GGH tables given in table 2 are included in the data card below) • 

$PROG KK=1 , FLIGHT=557, RUN=1 , NDZH«32, NGGH=24, 

DZHTABL=2 300 • , 175. ,5000. , 202. ,7000. , 240. ,9000. , 287. , 11000. , 

340. . 1 5000. . 450. . 20000. . 650. . 25000. . 772. . 3000. . 895. . 31000. , 

91 5. . 35000. .915.. 38000. . 900. .40000. . 850. . 42000. . 866. . 44000, 

,896. ,46000. ,920., 

GGHTABL=5000. ,0. ,0. , 1 1 000. ,0.5,30., 20000 . , 0. 9 , 27 . , 25000. , 

1 .1,35., 30000. ,1.25,45., 35000. ,1.2,40., 40000. ,1 .4,45., 

46000. . 1.9. 45., $ 

6/ 7/8/9 card 


Program Output 

A sample page of output is shown in figure 3. For this particular case, the 
level acceleration-deceleration method with the nonsurvey option was used. For 
each time point, angle of attack, angle of sideslip, indicated Mach number, indi- 
cated pressure altitude, stagnation pressure, and total temperature are listed, in 
addition to the corrections to Mach number, altitude, and static pressure due to 

position error. The corrections are labeled DM , DHP , and DPR , respectively. 

The last letter of each parameter specifies the method used to calculate the posi- 
tion error. If QQ = 0, appropriate headings for CP output are made. Any other 
integer for QQ results in appropriate headings for DHP output. 
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Conclusions 


The computer program presented calculates static pressure position error of 
airspeed systems. The program contains five separate methods for determining posi- 
tion error, of which the user may select from one to all five at a time. All the 
methods require an input data file consisting of aircraft and radar parameters. In 
addition, most of the methods require rawindsonde data or an atmospheric analysis, 
or both. 


National Aeronautics and Space Administration 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California, Jan . 9, 1984 


APPENDIX A — SUPERSONIC MACH NUMBER EQUATION 


The following derivation for calculating supersonic Mach number (PT/P > 1.893) 
was developed by Raymond Jackson of NASA Ames Research Center, Dryden Flight 
Research Facility. 


Let 


- J-iJL „2 

P 2 


r (i + y) 2 m 2 v(y-i) 

[4YM2 - 2(Y - 1). 


(eq. 3-26, ref. 1) 


1 + 


QCIT 

P 


1 + Y 
2 


M 2 


Y + 1 M 2 
2 


2Y 

_ Y + 1 


M 2 



Y z 1 

Y + 1 


) 


1/(Y-1) 


K 




and substitute Y = 1.4: 


K 


MC 2 (f ) 7 


0.18393711 


Let 


Y ~ 1 
2YM 2 


1 

7M 2 


Then 


1 


K 

QCIT 

P 


X(1 


x)1/(Y-1) = x ( 1 


x) 5/2 


Binomial expansion yields 


y 


K 

1 + QC IT 

P 



5_ 

16 


X 4 



3 v 6 _ 5 v 7 

256 X 1004 *•* 


and reversion of series gives 


x = y + 2.5y 2 + 10.625y 3 + 55y4 + 316.05469y5 
+ 1938y6 + 12421. 91 9y7 
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Let 


z 


lOy 


0,1839371 \ = 1 ,83971 
1 + QCIT ) " QCIT 

p / p 


Then 


x = 0. 1 z + 0.025z2 + 0.010625Z 3 + 0.0055z 4 

+ 0.0031605468Z 5 + 0.001938z 6 + 0. 001 2421 91 9z 7 

— = — - 2.5 - 0.4375Z - 0.175z2 - 0.08832031 5z 3 
X z 

- 0.05005z 4 - 0.030420265Z 5 
9 1 1.4285714 

M 2 = = 0.35714286 - 0.0625z 

- 0.025z2 - 0.01 261 71 87z 3 - 0.007l5z 4 

- 0.004345752z5 - 0.0087725z8 

The last term of this equation was added empirically so that the equation is exact 
at Mach 1. The effect of the empirical term is null at high Mach numbers. 

M = [(1.428572 - 0.3571428z - 0.0625z 2 - 0.025z 3 

- 0.01 2161 72z 4 - 0.00715z5 - 0.0043458z<5 

- 0.0087725z 9 )/z] 1/2 


where z = 1 . 839371 (P/PT) . 
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APPENDIX B - PERROR LISTING 


The following is a complete listing of the PERROR program. It is written in 
FORTRAN 77 and is run on an ELXSI computer at NASA Ames-Dryden. 



PROGRAM PERROR 


C 

C THIS PROGRAM CALCULATES POSITION ERRORS USING THE OPTIONS 
C OF SEVERAL DIFFERENT METHODS 

C 

DIMENSION DZHTABL(40) ,GGHTABL(60) 

DIMENSION EMINMAX(10),EAPffi0X(10),ZHPDIF(10), 

- ISTSV(4), IETSV(4), ISTAD(4), IETAD(4) 

INTEGER FLIGHT, RUN, 00, QQ 
DIMENSION IT(4) 

DATA LINE /6/ 

NAMELIST /PROG/ II, KK, LL, KW, NN, 00, QQ, IUNTS, HPREF, 

♦FLIGHT, RUN, DZ, NDZH, DZHTABL , NGGH, GGHTABL, ISURVEY, 

♦ISTSV, IETSV, l ST AD, IETAD 

DATA II, KK, LL, M, NN, 00, QQ, IUNTS, HPREF, FLIGHT, RUN, 

♦DZ, ISURVEY /I 4*0/ 

LINE=0 

OPEN (UNIT=1,FILE='perrorin') 

OPEN (UN I T=3 , F I LE= 'perrorout ' ) 

OPEN (UN IT=4,FI LE= ' f rmerge ' , FORM* ' UNFORMATTED ' ) 

C 

C READ NAMELIST CARDS, SOME OF WHICH ARE OPTIONAL DEPENDING ON METHODS USED 

C 

READ(1,PR0G) 

C 

C IF A SURVEY RUN IS BEING PROCESSED, READ FLIDAB-RADAR DATA DURING SURVEY, 
C CREATE DZEN TABLE, AND CALCULATE DZES. 

IF( ISURVEY .EQ .0) GO TO 7 

MSTSV * 3800000* I $TSV(1) + 60000* ISTSV(2) + 1000*ISTSV(3) + 

- I$TSV(4) 

METSV * 3600000* I ETSV(1) + 60000* IETSV (2) + 1000*IETSV(3) + 

- IETSV(4) 

MSTAD * 3600000* ISTAD(I) + 60000*ISTAD(2) + 1000*i$TAD(3) + 

- ISTAD(4) 

METAD = 3600000* I ETAD(1) + 60000*IETAD(2) + 1000*IETAD(3) + 

- IETAD(4) 

EMINMAX(I) = 1.0E10 
EMINMAX(IO) = -1.0E10 

30 READ(4,END=31) TIIyE,Z,R,E ,T,RRA,PR, 

- VCT,HPT,A,B,PT,PICT,QCIT,TT,FPA,FPH 
£*E/57.3 

RRA«RRA/57.3 

X=R*CQS(E)*$IN(RRA) 

Y*R*CO$(E)*CO$(RRA) 

IDTTL=TIME*1000. 

IF( IDTTL.LT. MSTSV) GO TO 30 
IFODTTL.GT. METSV) GO TO 31 
IF(E.LT.EMIW4AX(1)) EMINMAX(1)=E 
IF(E.GT.EMINMAX(10)) EMIMi4AX(10)*E 
GO TO 30 

31 DO 32 1*2,9 

EMINMAX(I) * EMINMAX(I) + (M)*((EMIMi(AX(10)-EMINMAX(1))/9.0) 

32 CONTINUE 
REWIND 4 

DO 33 1*1,10 
EAPPROX(I) - 1.0E10 
ZHPDIF(I) = 0.0 

33 CONTINUE 

34 READ(4,END*36) TIME,Z,R,E,T,RRA,PR, 

- VCT,HPT,A,B,PT,PICT,QCIT,TT,FPA,FPH 
E*E/57.3 
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RRA=RRA/57.3 

X=R*COS(E)*SIN(RRA) 

Y=R*CQS(E)*COS(RRA) 

IDTTL=TIME*1000. 

IF( IDTTL .LT.MSTSV) GO TO 34 
IF( IDTTL.GT.METSV) GO TO 36 
DO 35 1=1,10 

IF(ABS(E“EMINMAX(l)).GT.ABS(EAPPROX(l)-EMINMAX(l») GO TO 35 
EAPPROX(l) = E 
ZHPDIF(I) = Z-HPT 

35 CONTINUE 
GO TO 34 

36 DZES = ZHPDIF(IO) 

REWIND 4 

DO 37 1=1,10 

WRITE(3,101) EMINMAX( l),EAPPROX(l),ZHPDIF( I) 

101 FORMAT (1H .3F12.5) 

37 CONTINUE 


Of$G:NAt 

OF POOR 


C CALL SUBROUTINE TO WRITE TOP OF PAGE HEADINGS 

C 

7 CALL HEAD( FLIGHT, RUN) 

C ***** READ FL I DAB AND RADAR DATA FROM MERGED DATA BASE. 
10 READ(4,END=99) T)fc€,Z,R,E ,T,RRA,PR, 

-VCT ,HPT , A,B ,PT t P ICT ,QC IT ,TT ,FPA ,FPH 
E-E/57. 3 
RRA=RRA/57.3 
X=R*CO$(E)*$IN(RRA) 

Y=R*COS ( E ) *COS (RRA) 

IDTTL=TIME*1000. 

IFOSURVEY.NE.O .AND. IDTTL.LT.MSTAD) GO TO 10 
I F ( I SURVEY . NE . 0 .AND. IDTTL .GT .MET AD) GO TO 99 
IF(KK.EQ.O .AW. LL.EQ.O) GO TO 6 
C ***** OBTAIN DZH FROM TABLE DZHTABL. 

IF(Z.LE.DZHTABLO)) DZH=DZHTABL(2) 
IF(Z.GE.DZHTABL(NDZH-1)) DZB-DZHTABL(NDZH) 
IF(Z.LE.DZHTABLO) .OR. Z.GE.DZHTABL(NDZH-I)) GO TO 3 
DO 1 l*3,NDZH,2 
IF(Z.LT.DZHTABLO)) GO TO 2 

1 CONTINUE 

2 RATIO = (Z-DZHTABL( l-2))/(DZHTABL( I )-DZHTABL( (-2)) 

DZH - DZHTABL(I-I) + RATI0*(DZHTABL(l+1)-DZHTABL(i-1» 
IF(KK.EQ.O) GO TO 6 

IFOSURVEY.NE.O) GO TO 8 
C ***** OBTAIN G AND GH FROM TABLE GGhTABL. 

3 IF(Z.LE.QGHTABLO)) G=GGWTABL(2) 

IF(Z.LE.GGHTABLO)) GH-GGHTABL(3) 

I F(Z.GE. GGHTABL (NGGH-2)) G-GGHTABL(NGGH-I) 
IF(Z.GE.GGHTABL(NGGH-2)) GH-GGHTABL(NGGH) 
IF(Z.LE.GGHTABLO) .OR. Z.GE. GGHTABL (NGGH-2)) THEN 
G*G/6080.0 
GO TO 6 
END IF 

DO 4 l*4 f NGGH,3 
IF(Z.LT.GGNTABL(I)) GO TO 5 

4 CONTINUE 

5 RATIO - (Z-GGHTABL( l-3))/(GGHTABL( I )-GGHTABL( 1-3)) 

G - GGHTABL (1-2) + RAT 1 0* (GGHTABL (1+1 )-QGHTABL ( I -2) ) 

G - G/6080.0 

GH = GGHTABL(I-I) + RAT 1 0*(GGHTABL ( I +2)-GQHTABL( I -1 ) ) 
GO TO 6 


PACE IS 

quality 


***** OBTAIN DZEN FROM TABLES EAPPROX AND ZHPOIF. 

8 I F ( E . LE . EAPPROX ( 1 ) ) DZEN=ZHPD I F( 1 ) 

IF(E.GE.EAPPR0X(10)) DZEN-ZHPOIF(IO) 

IF(E.LE.EAPPROXd) .OR. E.GE.EAPPROX(IO)) GO TO 6 
DO 9 1=2,10 

IF(E.LT.EAPPROXd)) GO TO 11 

9 CONTINUE 

11 RATIO = (E-EAPPROX( 1-1 ))/(EAPPROX( I )-EAPPROX( 1-1 )) 

DZEN = ZHPOIF(M) + RAT IO*(ZHPD IF( l)-ZHPDIF(l-l)) 

TEST LINE COUNT AND CALL HEAD IF NECESSARY 

6 LINE = LINE + 1 

IF (M0D(LINE,6) .EQ. 0) CALL HEAD( FLIGHT, RUN) 

CALL SUBROUTINE UNITS TO MAKE UNIT CONVERSIONS TO ENGLISH UNITS IF NECESSARY 

IF( IUNTS.EQ.O) GO TO 15 

CALL UNITS ( IUNTS,PT,PICT,QCIT,TT) 

CALCULATE PT IF QCIT AND PICT OR VCT Ah© HPT ARE USED AS INPUTS 

15 IF(OO.EQ.O) GO TO 20 
IF(00.EQ.2) GO TO 17 

16 PT-QCIT+PICT 
GO TO 20 

17 CALL QCC(VCT,QCIT) 

CALL PRES(HPT,PICT) 

GO TO 16 

CALCULATIONS OF INDICATED MACH NO Ah© INDICATED PRESSURE ALTITUDE 

20 IF(PT.LT.PICT) GO TO 10 
PRATI-PT/PICT 

IF(OO.NE.O.AND.QQ.EQ.O) GO TO 21 
QCIT-PT-PICT 

21 CALL MACH(PRATI.AMI) 

CALL HPC (PICT, HP!) 

WRITE HEADING FOR FLIGHT COCITIONS, WRITE CONDITIONS AND HEADING FOR 
POSITION ERRORS 

WRITE (3,901) 

CALL SECTOT (TIME, IT) 

WRITE(3,902) IT, A, B, AMI ,HPI ,PT,TT 
IF(QQ.NE.O) GO TO 23 
WRITE (3,903) 

GO TO 24 
23 WRITE (3,916) 

POSITION ERROR CALCULATION USING RADAR-RAW I NSOCE METHOD 


24 IF ( 1 1 .EQ .0) GO TO 40 
P-PR 

PRAT-PT/PR 

CALL MACH (PRAT, AM) 

CALL PECALC(P,PICT,KIT,HP,H > I , AM, AMI ,DM,DPR,Dtf\CP,QQ) 
IF(QQ.NE.O) GO TO 25 
WRITE (3,904)DM,DPR,DH 5 
GO TO 40 

25 WRITE (3,910) DM,DPR,CP 


16 


c 

C POSITION ERROR CALCULATION USING LEVEL DECEL METHOD 

C 

40 IF(KK.EQ.O) GO TO 50 

I F( I SURVEY . EO . 0) CALL LEVDNS(AM,G,GH,HP,P,PT,Z,DZH,X,Y,RRA,DZ) 

IF( I SURVEY . NE . 0) CALL LEVDS(AM,HP,P,PT,Z,DZH,DZ,DZEN,DZES) 

CALL PECALC(P,PICT,QCIT,HP,HPI ,AM,AMI ,DM,DPR,Dtf\CP,QQ) 

IF(QQ.NE.O) GO TO 41 
WRITE(3,906)DM,DPR,DHP 
GO TO 50 

41 WRITE (3,912) DM, DPR, CP 

C 

C POSITION ERROR CALCULATION USING DESCENT PRESSURE METHOO 

C 

4 50 IF(LL.EQ.O) GO TO 60 

CALL DESCP(DZ,Z,DZH,P,PT,AM,HP) 

CALL PECALC(P,PICT,QCIT,HP,HPI ,AM,AMI ,DM,DPR,DHP,CP,QQ) 

IF(OQ.NE.O) GO TO 51 
WRITE(3,907) DM,DPR,DHP 
GO TO 60 

51 WR1TE(3,913) DM,DPR,CP 

C 

C POSITION ERROR CALCULATION USING DESCENT TE^ERATURE METHOD 

C 

60 IF(W.EQ.O) GO TO 70 

CALL DESCT(Z,HPREF,AMI ,TT,PT,P$,Ai£,HPS) 

CALL PECALC(PS,PICT,QCIT,HPS,HP1 , AMS, AM I ,DM,DPR,DHP,CP,QG) 
IF(QQ.NE.O) GO TO 61 
WRITE (3,908)DM,DPR,DHP 
GO TO 70 

61 WRITE (3,914) DM.DPR.CP 

C 

C POSIT IT ION ERROR CALCULATION USING TOTAL TEMPERATURE METHOD 

C 

70 IF (NN.EQ.O) GO TO 10 
CALL TOT(TT ,T,P,PT,AM) 

CALL PECALC(P,PICT,XIT,HP,HPI ,AM,AMI ,DM,DPR,DHP,CP,QQ) 

IF(QQ.NE.O) GO TO 71 
WRITE (3,909)DM,DPR,DHP 
GO TO 98 

71 WRITE (3,915) DM,DPR,CP 

98 GO TO 10 

99 CLQSE(I) 

CL0$E(3) 

CLQ$E(4) 

STOP 

C 

C FORMAT STATEMENTS 
C 

901 F0RMAT(1X > /,T3/TI1€(HR f MIN,SEC)',T25/ALPHA f DEG' f T45/BETA, DEG', 
•T63/IW MACH N0',T79,'I*C PRESSURE ALT' ,T99, 'STAG PRESSURE' ,T121, 
•'TOT TEW>'/T85, 'FEET', T104,'PSFA',T122, 'DEG, R'/) 

902 F0RMAT(1X, 17,3(12, IX), I3,T28,F5.1,T47,F5.1, 
*T67,F5.3,T85,F6.0,T103,F6.1,T122,F6.1//) 

903 FORMAT (1X,T12.'Wr ,T18, '«W',T25,'D^',T32,'DMLD',T38,'DPRLD', 

- T45, 'DHPLD' ,T52, 'OKP' J58, 'DPRCP' ,T65, 'DOTP' ,T72, 'HOT' , 

- T78, 'DPRDT' ,T85, 'DHPDT' ,T92, 'DMTT' ,T98, 'DPRTT' ,T105, 'DffTT'/) 

904 FORMAT (1X,T10,F6.4,T16,F6.3,T23,F7.0) 

906 FORMAT (V,T31,F6.4,T37,F6.3,T44,F7.0) 

907 FORMAT ( V,T51,F6.4,T57,F6.3.T64.F7.0) 
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908 FORMAT ( V ,T71,F6.4,T77,F6.3,T84,F7.0) 

909 FORMAT (V ,T91,F6.4,T97,F6.3J104,F7.0) 

910 FORMAT (1X,T1Q,F8.4,T17,F6.3,T24,F6.4) 

912 FORMAT (V ,T30,F6.4,T37,F6.3 ,T44,F6.4) 

913 FORMAT (V ,T50,F6.4,T57,F6.3,T64,F6.4) 

914 FORMAT (V ,T70,F6.4,T79,F6.3 ,T84,F6.4) 

915 FORMAT (V ,T90,F6.4,T97,F6.3,T104,F6.4) 

916 FORMAT (1X.T12, 'DWTJ18, 'DPRR' ,T26, 'CPR' ,T32, 'DM.D' ,T38, 'DPRLD' , 

- T46, 'CPLD' ,T52, 'DfcCP' ,T58, 'DPRDP' ,T66, 'CPDP' ,T72, 'DfcGT', 

- T78/DPRDT\T86/OTT'J92,'0Mn / ,T98/DPRTT' J106/CPTTV) 

END 

C 

SUBROUTINE SECT0T(t lme,lt) 

C 

C THIS SUBROUTINE CONVERTS A SINGLE REAL TIME WORD REPRESENTING 

C TOTAL SECONDS TO A 4 WORD INTEGER ARRAY REPRESENTING 

C hr,fflln,sec,msec 

C 

real mln.iBsec 
Integer it(3) 

C 

hr-tlme/3600. 

It(1)=hr 

C 

min-t lme/60. 

It(2)-mlrvlt(1)*60 

C 

seotlme 

it(3)=sec- lt(2)*60- 1 1(1 )*3600 
C 

msec-tlne* 1000. 

It(4)««sec-lt(3)*1000-lt(2)*6.e04-lt(1)*3.6e06 

C 

return 

end 

C 

C 

SUBROUTINE HE AD (FLIGHT, RUN) 

C 

C THIS SUBROUTINE INRITES MAIN HEADINGS AT THE TOP OF EACH PAGE 

C 

NR ITE(3, 100)FL IGHT ,RUN 

100 F0RMAT(1H1,T36, 'STATIC PRESSURE POSITION ERROR CORRECTIONS FOR FLI 
*GHT M3,' RUN ',12,///) 

RETURN 

END 

C 

SUBROUTINE UNITS(IUNTS,PT,PICT,QCIT,n) 

THIS SUBROUTINE CONVERTS PARAMETERS IN METRIC UNITS TO ENGLISH UNITS 

IF0UNTS.EQ.2) GO TO 5 
PT*144.*PT 
PICT»144.*PICT 
QCIT«144.*QCIT 
GO TO 10 
5 PT-.020885*PT 
PICT«.020885*PICT 
QC I T». 020885*00 IT 
TT - 1 . 8*(TT-459 . 67)+273 . 1 6 
10 RETURN 
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OF POOR QUALITY 

ENO 

c 

SUBROUTINE PECALC (P.PICT.QCIT.HP.HPI , AM, AM I .OM.DPR.OHP.CP.QQ) 

C 

C THIS SUBROUTINE CALCULATES POSITION ERRORS COMMON TO ALL OF THE 
C SUBROUTINES 

DPR- (P— P ICT)/P 
DM- AM- AM I 

IF (QQ.NE.O) GO TO 5 
CALL HPC(P.HP) 

DHP-HP-HPI 
GO TO 10 

5 CP*(P-P!CT)/QCIT 
10 RETURN 
END 

C 

SUBROUTINE MACH (PRAT, AM) 

C 

C THIS SUBROUTINE CALCULATES MACH NUlfcER AS A FUNCTION OF PRESSURE RATIO 

C 

IF (PRAT.GT. 1.89293) GO TO 5 
AM-SQRT(5*((PRAT)*V2857-U) 

GO TO 10 

5 Y-1. 839371/PRAT 

ASSORT ( ( 1 . 42857- . 3571 43*Y- . 0625*( Y**2)- . 025*( Y**3)- . 01 261 7*( Y**4) 
*-.00715*(Y**5)-.QQ43458*(Y**6)-.QQ87725*(Y**9))/Y) 

10 RETURN 
END 

C 

SUBROUTINE PRAT 10 (AM.PRAT) 

C 

C THIS SUBROUTINE CALCULATES PRESSURE RATIO AS A FUNCTION OF MACH NUMBER 

C 

IF(AM.LT.1.) GO TO 5 

PRAT-1 . 2*AM**2*( ( (5 . 76*AM**2)/( 5 . 6*AM**2- . 8))**2 . 5) 

GO TO 10 

5 PRAT-(1.+.2*AM**2)**3.5 
10 RETURN 
END 

C 

SUBROUTINE HW.HP) 

C 

C THIS SUBROURINE CALCULATES PRESSURE ALTITUDE AS A FUNCTION OF STATIC 
C PRESSURE 

IF(P.LT.114.35) GO TO 5 

IF(P.LT.472.68) GO TO 10 

HP— 145442. *( (P/21 16.22)**. 190262-1 . ) 

GO TO 15 

5 HP-6561 6. 8*71 0794. *(( . 008746*P)**(-. 029271 )-1 .) 

GO TO 15 

10 FP-36089 . 2-20805 . 7*( ALOG( . 0021 1 55*P ) ) 

GO TO 15 
15 RETURN 
END 

C 

SUBROUTINE PRES (tt>,PS) 

C 

C THIS SUBROUTINE CALCULATES STATIC PRESSURE AS A FUNCTION OF PRESSURE 
C ALTITUDE 

I F (HP. LT. 38089. 2) GO TO 5 
IFOP.LT.65616.8) GO TO 10 



P$=51 . 9754* (( 1 - . 0000042204* (HP-6561 6 . 8) )**(-!> . 38/ttv , 

GO TO 15 

5 PS*2116.22*(1 .-.0000068754*tf > )**5.2561 
GO TO 15 

10 PS«472. 68*(2 . 71828**(- . 000048063* (HP-36089 . 24) ) ) 

15 RETURN 
END 

C 

SUBROUTINE QCC(VCT,QC) 

C 

C THIS SUBROUTINE CALCULATES QC AS A FUNCTION OF VCT 

C 

IF(VCT.GT.661.4) GO TO 5 

X*21 16. 22* ((.1 42857* (. 0000032*VCT**2)+1 .) **3.5-1) 

GO TO 10 

5 QC* . 005805* VCT**2*(5 . 76/(5 . 6-349964 . /VCT**2) )**2 . 5-21 1 6 . 22 
10 RETURN 
END 

C 

SUBROUTINE LEVONS (AM.G.GH.HP.P.PT.Z^OZH.XJ^A.DZ) 

C 

C 

C THIS SUBROUTINE CALCULATES TRUE MACH NUMBER AND STATIC PRESSURE 
C BY THE LEVEL DECELERATION METHOD FOR NON-SURVEY TEST. 

C 

DANG * ABS(RRA-GH)/57.3 
DR * SGRT(X**2 + Y**2) 

DHPG = DR*G*COS(DANG) 

HP - Z-DZH+OHPG-OZ 
CALL PRES(tf\P) 

PRAT-PT/P 
CALL MACH(PRAT.AM) 

RETURN 

EM) 

SUBROUTINE LEVDS(AM,HP,P,PT,Z,DZH,DZ,DZEN,DZES) 

C THIS SUBROUTINE CALCULATES TRUE MACH NUMBER AM) STATIC PRESSURE 

C BY THE LEVEL DECELERATION METHOO FOR SURVEY TEST. 

HP - Z-DZH-DZEN+DZES-DZ 
CALL PRES (HP»P) 

PRAT - PT/P 
CALL MACH (PRAT, AM) 

RETURN 

END 

SUBROUTINE DE$CP(DZ,Z,DZH,P,PT,AM,HP) 

C THIS SUBROUTINE CALCULATES TRUE MACH NUMBER AND STATIC PRESSURE 
C BY THE DESCENT STATIC PRESSURE METHOD 
C 

HP - Z-OZ-DZH 
CALL PRES(HP.P) 

PRAT - PT/P 
CALL MACH(PRAT,AM) 

RETURN 

EM) 

SUBROUTINE TOT(TT ,T,P,PT,AM) 

C 

C THIS SUBROUTINE CALCULATES MACH NUMBER FROM kCASUREMENTS OF TOTAL 
C TEMPERATURE AM) AMBIENT TEMPERATURE 
AM-SQRT(5*(TT-T)/T) 
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original page ; IS 

CALL PRAT I0(AM, PRAT) OF POOR QUALITY 

P-PT/PRAT 

RETURN 

END 

SUBROUTINE DESCT(Z,HPREF,AMI ,TT,PT,P,AM,HP) 

C THIS SUBROUTINE CALCULATES TRUE MACH NUMBER AND STATIC PRESSURE 
C BY THE DESCENT TEMPERATURE METHOD 

C 

C STATEMENT FUNCTION FOR CALCULATING AMBIENT TEMPERATURE 

C 

LOGICAL FIRST 
DATA FIRST/. TRUE./ 

IF(. NOT. FIRST) GO TO 10 
FIRST - .FALSE. 

HP4PREF 
CALL PRES(HP.P) 

PRAT-PT/P 
CALL MACH(PRAT,AM) 

GO TO 60 
C 

C CALCULATION OF TRUE PRESSURE ALTITUDE USING ITERATION 

C 

10 DZ-Z-ZJ 

c 

C SET COUNTER FOR NUMBER OF ITERATIONS LIMITED TO 10 

C 

N-0 

HP04P 

ZBAR « (Z+ZJ)/2. 

TTBAR-(TT+TTJ)/2. 

AMD- AM I 

DH-(20825000.*DZ)/(20825000.+ZBAR) 

C CALCULATION OF AMBIENT TEMPERATURE AS A FUNCTION OF ALTITUDE 

C 

C 

IF(ZBAR.GT. 36089.) GO TO 20 
TS-518.7*(1 ,-.0000068756*ZBAR) 

GO TO 40 

20 IF(ZBAR.GT. 65617.) GO TO 30 
TS-390. 

GO TO 40 

30 TS-390. *(!.+. 000001 4069*(ZBAR-6561 7.)) 

40 TA-TTBAR/(1 .0 + 0.2*AMC**2) 

DFP-TS*DH/TA 
HP4P0+0MP 
CALL PRES(HP,P) 

PRAT-PT/P 
CALL MACH(PRAT.AM) 

DM I -AMI- AMD 

IF(ABS(DMI).LT.0.0005) GO TO 60 

AMC-AM 

N-Kfl 

IF(N.GT.IO) GO TO 60 
GO TO 40 
60ZJ-Z 
TTJ-TT 
RETURN 
END 
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TABLE 1 . — NECESSARY METHOD INPUTS TO PERROR 


LEVD 



Survey 

( ISURVEY = 1 ) 

Nonsurvey 
(ISURVEY = 0) 

DESCP 

DESCT 

TOT 

RADAR 

$ PROG card 

KK = 1 

KK = 1 

LL = 1 

MM = 1 

NN = 1 

II = 


ISTSV, IETSV 

DZHTABL 

DZH 

HPREF 




ISTAD , IETAD 

DZ 

DZ 





DZHTABL 

NDZH 

NDZH 





DZ 

NGGH 






NDZH 






Merged file 

P 

P, PT* 

P 

P 

P 

PR 


PT* 

GGHTABL 

PT* 

PT* 

PT* 

PT* 


Z 

Z, RRX, RRY 

Z 

Z 

TT 

Z 



RRA 


TT 



Calculated 

DZEN 

DANG 

MI 

MI 

MI 

MI 


DZES 

DR 



T 



MI 

MI 






*pCIT and PICT or HPT and VCT may be substituted. 


TABLE 2. — ATMOSPHERIC ANALYSIS RESULTS 
FOR F-1 4 FLIGHT 557 


Z, 103 ft DZH * 

G, ft/nmi 

GH, deg 
true north 

46 

920 

1.9 

045 

44 

896 

— 

— 

42 

866 

— 

— 

40 

850 

1 .4 

045 

38 

900 

— 

— 

35 

915 

1.2 

040 

31 

915 

— 

— 

25 

772 

1.1 

035 

20 

650 

0.9 

027 

15 

450 

— 

— 

1 1 

340 

0.5 

030 

9 

287 

— 

— 

7 

240 

— 

— 

5 

202 

0.0 

— 

surface 

175 

— 

— 
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Fllghtpath 


RRX, RRY 

(X, Y) 



(a) Nonsurvey method . 



(b) Survey method . 

Figure I. Methods used to compensate 
for the horizontal pressure gradient 
when using level acceleration- 
deceleration . 



Figure 2. Block diagram of 
input sources to PERROR • 
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STATIC PRESSURE POSITION ERROR CORRECTIONS FOR FLIGHT 178 RUN 1 


TIME(HR,MIN,SEC) 

ALPHA, DEG 

BETA, DEG 

DC MACH NO 

110 PRESSURE ALT 
FEET 

STAG PRESSURE 
PSFA 

TOT TEkP 
DEG. R 

8 23 30 0 

6.7 

***** 

.780 

27851. 

1035.3 

-7.1 

OMR DPRP 

DtffR DMLD DPRLD 

DHPLD DM3P 

DPRDP DHPDP DM)T 

DPRDT ttfDT DMTT 

OPRTT (H>n 



.0197 -.020 

433. 





TIME (lit, MIN, SEC) 

ALPHA, DEG 

BETA, DEG 

1 tO MACH NO 

ItO PRESSURE ALT 
FEET 

STAG PRESSURE 
PSFA 

TOT TEM> 
DEG, R 

8 23 30 500 

6.7 

***** 

.780 

27851. 

1035.3 

-7.1 

DM1 DPRP 

DHPR DMLD DPRLD 

DHPLD DMDP 

DPRDP OtPOP DMDT 

DPRDT DHPDT DMTT 

DPRTT DMTT 



.0196 -.019 

430. 





TIME(Ht,MIN,SEC) 

ALPHA, DEG 

BETA, DEG 

It© MACH NO 

ItO PRESSURE ALT 
FEET 

STAG PRESSURE 
PSFA 

TOT TEM> 
DEG, R 

8 23 31 0 

6.8 

•••** 

.783 

27862. 

1037.6 

-6.9 

DMR OPRP 

DtfR [HO DPRLD 

OtfL D OKP 

DPRDP DHPOP DM)T 

DPRDT DHPDT DMR 

DPRR DHPR 



.0193 -.019 

426. 





T 1 ME ( HR , M 1 N , SEC ) 

ALPHA, DEG 

BETA, DEG 

IND MACH NO 

ItO PRESSURE ALT 
FEET 

STAG PRESSURE 
PSFA 

TOT TEM> 
DEG, R 

8 23 31 500 

6.8 

***** 

.783 

27862. 

1037.6 

-6.9 

DM1 DPR P 

DHPR MS) DPRLD 

DHPLD DM? 

DPRDP OHPDP DM)T 

DPRDT DHPDT DMn 

DPRR DHPR 



.0196 -.020 

431. 





TIME (HR, MIN, SEC) 

ALPHA, DEG 

BETA, DEG 

ItO MACH NO 

ItO PRESSURE ALT 
FEET 

STAG PRESSURE 
PSFA 

TOT TEM> 
DEG. R 

8 23 32 0 

6.6 

***** 

.786 

27873. 

1039.6 

-6.7 

DM1 DPRP 

DHPR DM LD DPRLD 

DfrPLD DM3P 

DPRDP DHPOP DM)T 

DPRDT DtPDT DMn 

DPRR DHPR 



.0194 -.019 

428. 






Figure 3. 

Output 

format of PERROR program . 
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